Generating "squeezed" superpositions of coherent states using photon addition and 

subtraction 
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We study how photon addition and subtraction can be used to generate squeezed superpositions of 
coherent states in free-traveling fields (SSCSs) with high fidelities and large amplitudes. It is shown 
that an arbitrary iV-photon subtraction results in the generation of a SSCS with nearly the perfect 
fidelity (F > 0.999) regardless of the number of photons subtracted. In this case, the amplitude of 
the SSCS increases as the number of the subtracted photons gets larger. For example, two-photon 
subtraction from a squeezed vacuum state of 6.1dB can generate a SSCS of a = 1.26, while in the 
case of the four-photon subtraction a SSCS of a larger amplitude a = 1.65 is obtained under the 
same condition. When a photon is subtracted from a squeezed vacuum state and another photon is 
added subsequently, a SSCS with a lower fidelity (F ~ 0.96) yet higher amplitude (a w 2) can be 
generated. We analyze some experimental imperfections including inefficiency of the detector used 
for the photon subtraction. 

PACS numbers: PACS number(s); 42.50.Dv, 03.67.-a, 42.50.Ex 



I. INTRODUCTION 

The development of the quantum theory of light has 
deepened our understanding of nonclassical properties of 
optical fields. Recently, superpositions of coherent states 
in free-traveling fields (SCSs) (HQ have attracted special 
attention due to their remarkable usefulness. When their 
amplitudes are large, the SCSs show typical properties 
of macroscopic quantum superpositions, and because of 
this, they are often called "Schrodinger cat states" recall- 
ing the famous cat paradox The SCSs enable one to 
perform many interesting studies for fundamental tests of 
quantum theory 4, 5, 6]. Furthermore, it has been found 
that SCSs are useful for various applications in quantum 
information processing @, H, [tj [HI [ill EH • The power of 
this approach lies in the fact that all the four Bell states 
can be discriminated in a deterministic way only using a 
beam splitter and photon counting [H, Q , which is obvi- 
ously not the case for the single-photon based approach. 

In spite of the manifold usefulness of the SCSs, until 
recently, the generation of free-travelling SCSs has been 
known to be difficult. There have been schemes to gener- 
ate such SCSs using strong nonlinear interactions JJ, EH 
or photon number resolving detectors [Til [l5| , which are 
not feasible using current technology. Recently, more re- 
alistic schemes have been suggested by several authors 
[U, EE El 01, HI EH- For example, a scheme using 
weak Kerr nonlinearities and simple optical elements was 
suggested [l8[ based on a previous proposal where strong 
Kerr nonlinearities are required [13j |. As another exam- 
ple, a simple observation was made that SCSs with small 
amplitudes, such as a < 1.2, are well approximated by 
squeezed single photons It was also pointed out 

that squeezed single photons can be obtained by sub- 
tracting or adding one photon from pure squeezed vacu- 



ums [22j. Meanwhile, single-photon-subtracted squeezed 
states, which are close to the SCSs with small amplitudes 
( a ~ l)i have be en g enerated by several experimental 
groups 
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|25|. |26H and theoretical analysis has been 
13 |2g|. Recently, squeezed SCSs (SS- 
CSs) were generated and detected [30|, where the size 
of the states (a = 1.6) was reasonably large for funda- 
mental tests of quantum theory and quantum informa- 
tion processing, for which the states are suited despite 
their squeezing [3l[ . A scheme using time separated two- 
photon subtraction was suggested [32| and experimen- 
tally demonstrated .33] to generate SCSs of large am- 
plitudes. Despite all the recent progress, however, the 
fidelities of the generated states are yet to be improved 
for practical quantum information processing. 

The directions of the development for the generation 
of SCSs are twofold. First, one needs to generate SCSs 
with larger amplitudes (a > 2) for macroscopic tests of 
quantum theory. Second, for quantum information pro- 
cessing, it is important to generate SCSs with higher fi- 
delity F > 0.99 while a « 1.6 is an appropriate value 
[34[ ■ The SSCSs, generated in a recent experiment [30l |. 
are simply a squeezed version of the SCSs. Interestingly, 
the direction of the squeezing in Ref. [3fJ makes the SS- 
CSs more robust against decoherence than the regular 
SCSs (35| . The SSCSs can be useful in some protocols as 
they are [3l|, and if required, it may be possible to un- 
squeeze them by means of the squeezing transformation 
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In this paper, we are interested in finding methods to 
generate SSCSs with larger amplitudes and higher fideli- 
ties using photon subtraction and addition. Here, we 
show that the two photon subtraction enables one to pro- 
duce the SSCS with very high fidelity, and the combina- 
tion of the subtraction and addition can return the SSCS 
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with a lower fidelity yet higher amplitude. 

We also find that consecutive applications of photon 
subtraction (or subtracting a well defined number of pho- 
tons) from a squeezed vacuum state result in the genera- 
tion of a SSCS with nearly the perfect fidelity regardless 
of the number of photons subtracted. The amplitude 
of the SSCS increases as the number of the subtracted 
photons gets larger. This paper is organized as follows. 
In Sec. II, we investigate combinations of the ideal single 
photon addition and subtraction. In Sec. Ill, some exper- 
imental imperfections are analyzed for the realization of 
the states discussed in Sec. II. In Sec. IV, we numerically 
show that an arbitrary N photon subtraction, regardless 
of N, can be used to generate a SSCS with an extremely 
high fidelity as F > 0.999. We conclude with final re- 
marks in Sec. V. 



II. PHOTON ADDITION TO AND 
SUBTRACTION FROM SQUEEZED VACUUM 

A. Combinations of photon addition and 
subtraction 

The SCSs are defined as 

\SCS V ) = N v {\a) +e l *\ -a)), (1) 

where J\f v is a normalization factor, | ± a) is a coherent 
state of amplitude ±a, and <p is a real local phase factor. 
The SCSs such as |SCS±(a)) = Af±(\a)±\-a)) are called 
even and odd SCSs respectively because the even (odd) 
SCS always contains an even (odd) number of photons. 
The squeezed vacuum state, |S(r)), with the squeezing 
parameter r can be obtained by applying the squeezing 
operator, 

S(r) = e i( a2 - at2 ), (2) 

where a (aJ) is the bosonic annihilation (creation) oper- 
ator, to the vacuum state. In the number state basis, it 
can be represented as 

|S(r)) = V^x&jr V )^ [~ ^tanhr] fc |2fc), (3) 

fc=0 

where r was assumed to be real. As seen in Eq. (|3"T). 
the squeezed vacuum state contains only even number 
of photons. A squeezed single photon, known as a good 
approximation of a small SCS, can be obtained by adding 
a photon to a squeezed vacuum as 

tfS{r)\0) = cosh rS(r) 1 1), (4) 

where the right hand side is unnormalized due to the 
characteristics of the creation operator a* . We note that 
Eq. ((4]) can easily be shown using following unitary trans- 
formations (37l ]: 

S'(r)aS(r) — a cosh r — a* sinhr, 
S' (r)a* S(r) = a' cosh r — a sinh r. 



It is also known that the squeezed single photon can also 
be obtained by subtracting a photon from a squeezed 
vacuum as 

aS(r)\0) = -sinh rS(r) | 1). (6) 

This may cause us to conjecture that when the photon 
addition and subtraction are applied successively to the 
squeezed vacuum, an approximate even SCS may be gen- 
erated. We can consider four immediate cases, namely, 
addition and subtraction (aw), subtraction and addition 
(a* a), successive additions ((a*) 2 ), and successive sub- 
tractions (a 2 ). From Eqs. (|4]) and ([6]), it is straigtforward 
to see that (c^) 2 will result in the same state produced 
using aJa. It is also easy to see that a& will cause the 
same effect with a 2 . Therefore, we shall consider only 
two cases, a' a and a 2 , among the four based on the fact 
that the photon subtraction is relatively easier to perform 
than the photon addition. 

Now suppose an ideal situation that a photon is sub- 
tracted from the squeezed vacuum state and then another 
photon is subsequently added. The resulting state, which 
we shall call photon-subtracted-and-added squeezed state 
(PSAS), is obtained by applying the annihilation and 
creation operators, tra, to the squeezed vacuum state 
jS'(r)). After a straightforward calculation using Eqs. ([5]) 
and the normalization, the PSAS appears to be 

l^ata) =A4ta£(r)(|0) - v^tanhr)" 1 ^))) (7) 

where A/gta = {1 + 2(tanhr) -2 } -1 / 2 . In the same man- 
ner, the two-photon subtracted squeezed state (TPSS) 
can be obtained as 

1^) = AA a2 5(r)(|0> - V2tanhr|2») (8) 

with Af' a 2 = { 1 + 2(tanhr) 2 }- 1 /2. 

B. Fidelities against ideal states 

The fidelity, F = | (ip\ipt) 1 2 , is a measure of how close a 
state \ip) is to the target state \i/>t)- It is unity when the 
two states are identical, while it is zero when the two are 
orthogonal to each other. The fidelity between the TPSS 
(or PSAS) and the ideal squeezed (or regular) SCS can 
be obtained as follows. A even SSCS can be expressed as 

|SSCS) =N+S(r')[\a) + \ -a)], (9) 

where A/+ = [2 + 2e~ 2Q2 ]" 1/2 . The fidelity F between 
the TPSS and the even SSCS can be calculated as 

F= I I SSCS) 1 2 

- 8|A/- + A/- a2 [(0|£t( r _ r ')\a) + v{2\&(r - r')\a)}\ 2 

(10) 

where v = — tanhr. This expression can be evaluated 
with help of the x-representation: 

/oo 
{iPx\x){x\^)dx, (11) 
-OO 
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FIG. 1: (a) Optimal fidelity F opt between the ideal TPSS (a 2 ) 
and the squeezed SCS of amplitude a (solid curve) and the 
optimal fidelity between the ideal TPSS and the correspond- 
ing regular SCS (dashed curve). The squeezing parameter 
r' opf of the initial squeezed state is r — —0.7. (b) The squeez- 
ing parameter r' of the target SSCS for which the fidelity is 
optimized. 



FIG. 2: (a) Optimal fidelity F opt between the ideal PSAS 
(a^a) and the SSCS of amplitude a (solid curve) and the op- 
timal fidelity between the ideal PSAS and the corresponding 
regular SCS (dashed curve). The squeezing parameter r' opt of 
the initial squeezed state is r — —0.7. (b) The squeezing pa- 
rameter of the target SSCS for which the fidelity is optimized. 



where the relevant wave functions are 

(x\S(r-r')\0) = (Vr 1/4 e~ x2/2s2 

2a; 2 



(x\S(r-r')\2) = (W)- 1/4 ( 
(x\a) 



- 1 



= 7r -l/4 e -( K -v^ Q ) 2 /2 



(12) 



where g = exp[— (r — r')]. After some calculations, we 
arrive at the fidelity: 
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Af+Af a 2{l + v 
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1 + 4a z g 



(1+5 2 ) 



2\2 



} ■ (13) 



The same approach can be used to derive fidelity of states 
prepared by combined photon subtraction and addition: 



where 



,r ,r r l + 4a 2 g 2 - o 4 -, 2 



(1 



[l + 2 M 2 ]-^. 



li - 



tanh i 



(14) 



(15) 



In Figs. Q] and [21 we have used the analytical expres- 
sions to show the optimal fidelities for a range of co- 
herent amplitudes. That is, for each a, such squeezing 
r' of the SSCS is found to provide the maximal over- 
lap with the state prepared from initial squeezed state 
with r = —0.7 (p-squeezed state, approximately 6.1 dB 
of squeezing, which can be realized using current technol- 
ogy). For comparison a fidelity with the corresponding 
regular SCS is also shown. 



Figure QJa) shows the fidelity of the TPSS. It is seen 
that although the fidelity tops at 0.9 for a regular SCS, 
when a SSCS is considered, the fidelity of F = 0.999 can 
be achieved for a = 1.26 and r' — —0.425. The squeez- 
ing parameter r' of the target SSCS that optimizes the 
fidelity against the corresponding amplitude a is plotted 
in Fig. 1(b). The fidelity for the PSAS and the optimiz- 
ing squeezing parameter against a are depicted in Fig. [21 
The optimal fidelity is F = 0.956, which is not as good as 
the case of the TPSS. However, in this case, the ampli- 
tude of the SSCS is larger as a — 1.93, while the optimal 
squeezing for the SSCS is only r' = —0.17. 

The Wigner function of a state with density opera- 
tor p can be obtained from the Fourier transform of its 
characteristic function C(£) = Tr[D(()p], where D(() = 
exp[(a^ — (*a\ is the displacement operator. The Wigner 
functions of the TPSS and PSAS can be obtained using 
Eqs. (0 and © as 



W#{p)=N tl 



n -2Z 



1 + 2tanhr{Z' 
+ tanhr[l + 8Z{Z - 1)]} , (16) 
W^a(fi) = N as e- 2Z [\ + 2cothr{Z' 

+ cothr[l + 8Z{Z- 1)]}1 (17) 

where Z = e 2r /3 2 + e- 2r (3f and Z' = -4e 2r ff + Ae~ 2r (3f. 
In Figs. [3] and [H we consider the squeezing parameter 
r = —0.7 (6.1dB), for the initial squeezed state. Fig- 
ure [3] shows again that the TPSS is an extremely good 
approximation of the even SSCS. In Fig. [3J the Wigner 
functions of the TPSS of r = -0.7 and the even SSCS 
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(a) 




(b) 




FIG. 3: (a) The Wigner function of the TPSS with r = -0.7 
(«6.1dB) and (b) the Wigner function of the ideal SSCS with 
a = 1.26 and r' = -0.425. The fidelity between the two 
states is nearly perfect as F > 0.999. 



of r' = —0.425 look virtually identical and the fidelity 
between the two states is F > 0.999. Figure [4] presents 
the Wigner functions of the PSAS of r = —0.7 and the 
even SSCS of r' = —0.14 and a = 2, where the fidelity be- 
tween the two states is F w 0.955. This shows that slight 
variation of parameters still allows for high fidelity. 



To conclude, if one is interested in generating SSCSs 
of high fidelity, the two-photon subtraction would be a 
useful scheme, while the photon subtraction and addition 
would be a better strategy to generate large SSCSs. 



III. EXPERIMENTAL CONSIDERATIONS 

A. Optical operations with ideal avalanche 
photodetectors 

So far, we have considered ideal photon addition and 
subtraction using annihilation and creation operators. In 
real experiments, however, they can be implemented only 
using approximate schemes. The setup for implementa- 
tion of the annihilation (creation) of a single photon con- 
sists of a beam splitter - BS (noncollincar optical para- 
metric amplification - NOPA) and an avalanche photode- 
tector. In the Wigner function formalism, the actions 
of BS and NOPA, coupling two modes of light, can be 
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FIG. 



4: (a) The Wigner function of the PSAS (left) with 
r = —0.7 (~6.1dB) and (b) the Wigner function of the ideal 
SSCS with a = 2 and r' = -0.14 (right). The fidelity between 
the two states us F m 0.955. 



characterized with help of a transformation matrix act- 
ing on a vector of variables {x±,pi,X2,P2) corresponding 
to quadrature operators Xj and Pj with [Xj,Pj<] = i8jj>: 



t r 

t x 

-r t 

o -% o t 



where the generalized parameters are 



t = VT 



= x = VT^r 



(18) 



(19) 



for V12 to describe the action of a beam splitter with 
transmissivity T, and 



t = VG, r = - X = VG^l 



(20) 



for V12 to describe a NOPA with amplification gain G. 
The three-mode Wigner function for the initial squeezed 
state and two vacuum ancillas is expressed as 



Wtat(0 



1 



(2^) 2 ^E^7 



exp 



" 2^tot 



i,t\ 



(21) 



Here, the vector of variables for the three-mode Wigner 
function is defined as £ = {x\,pi, XAiPA, %b-,Pb)-> where 
x = (3 r /\/2 and p = /3i/v2 being compared to the vari- 
able, z, used in Eqs. (fTB)) and (fT7|) . and the subscripts of 
x and p in order to denote the initial and two ancillary 
modes, respectively. Furthermore, |.| denotes the deter- 
minant of the matrix, T stands for transposition and E to t 
is the covariance matrix of the state, 

E tot =diag(K i y p ,i i ii), (22) 
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with V x and V p being the variances of the initial squeezed 
state. 

The initial state subsequently interacts with the two 
vacuum modes, transforming the vector of variables into 



where 

Wx(x,p) 



2- 



exp 



1 



-(x,p)T, 1 {x,p) 1 



(31) 



= V B V A Z, 



(23) 



where the transformation matricies, Va and Vb are of 
the form (fT8|) coupling modes 1 and A and 1 and B with 
parameters tA, ta, XA and ts, r B and xb, respectiv- 
elly. To complete the transformation, a conditionning 
measurement is needed. When using two ideal avalanche 
photodetectors and post-selecting the state only when 
both produce a detection event, we implement a pair of 
projection operators 1 — |0)(0| and the Wigner function 
of the output state is transformed to: 

W out (x 1 ,p 1 ) = J W tot (£)[l-2TrW vac (x A ,PA)] 

x [1 - 2TrW mc (x B ,PB)]dxAdpAdx B dpB 

where W va . c (x,p) = exp(— x 1 —p 2 )/ir is the Wigner func- 
tion of a vacuum state. This integral can be expressed 
as a sum of four Gaussian integrals, 

1 



is the Wigner function of a Gaussian state with covari- 
ance matrix E, the coefficients are 



(32) 



and 



M = (c, 



2c, 



ny 



2c yn + 4c TO )(V|Stot|)" 1 (33) 



is the normalization factor and overall probability of suc- 
cess. 

In order to compare the final state (|30[) with our tar- 
get state, the even SSCS, we need to employ its Wigner 
function 



W SC s{x,p) = Ns GS - 
(x + \/2ga) 



-{exp 



[x - jTgaf 



7(S n „) ~ 2J(S nl( ) 



2ir^/\^\ 

- 2I{H yn ) + 4J(E ra ) 
where 7(E) is a shorthand for 



(25) 



exp • = - +2e~ x /9 cos(2y/2gap)}, 

(34) 

where Afscs = [2 + 2e~ 2 ""]~ 1 / 2 , a is the coherent am- 
plitude, and g = exp(— 2r) characterizes squeezing of the 
state. The fidelity can then be calculated as the overlap 
between the two Wigner functions: 



7(E) = 



2*) 



1 



exp ( ~2^ S £ I dxAdpAdx B dp B 



F = 2ir / W out (x,p)W S cs(x,p)dxdp. 



and Sjjv is a covariance matrix for the particular event 
when either the two detectors did detect vacuum (y) or 
no measurement has taken place (n). The particular co- 
variance matrices can be found as 



(35) 



Since the Wigner function W ou t(x,p) is a sum of Gaus- 
sian functions, we can treat the integration in parts and 
express the final fidelity as 



F 



1 



^jj' — Ptot "I" 11 jj'J 



n, 



(27) 



where the Hjji is the semi-inverted covariance matrix of 
the vacuum state, 

II„„ =diag(0, 0,0, 0,0,0), U yy = diag(0,0,2,2,2,2), 
U yn =diag(0,0,2,2,0,0), U ny = diag(0, 0, 0, 0, 2, 2). 

(28) 

The matrices can be decomposed into a block form, /*[diag(A*j v )\ —2 



MV^\ [Cnnf(Ann) 2Cnyf(Any) (36) 

_ 2Cy n f (A yn ) + 4:Cyyf(Ayy)^ . 

Here /(E) denotes the fidelity between a Gaussian state 
with covariance matrix E and a SSCS. Since all the co- 
variance matrices used are diagonal, we can write these 
partial fidelities as functions of the diagonal elements as 



(1 



1 



)J{g + 2p){- + 2v) 
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, , (29) 

'33' J W 

where the submatrix Ajji corresponds to variables Xi,pi, 
submatrix Bjji to xa,Pa, xb,Pb and submatrix Cjji cov- 
ers relations between these two groups. After the inte- 
gration, we can arrive at the final Wigner function: 

1 



exp( 



-2a 2 



1 + 2/V-? 



exp(- 



2a 2 



1 + l/2gv J 



W out (x,p) 



c nn W Ann {x,p) - 2c ny W Any (x,p) 



AA^E 

- 2c yn W Ayrl (x,p) + Ac yy W Ayy (x,p) 



(30) 



(37) 

With Eqs. (|3T>)) and (f3"T)) we can finally obtain the re- 
quired fidelity. In analogy with the previous section, the 
behavior of the state is depicted on Figs. [5] and [6] The 
comparison reveals that these real states behave in a simi- 
lar pattern as the ideal ones, the performance is however 
slightly worse. For realistic parameters, T = 0.99 and 
G = 1.01, the optimal fidelities are F w 0.96 for two 
photon subtraction and F sa 0.91 for subtraction and 
addition. 
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FIG. 5: (a) Optimal fidelity F opt between the TPSS (a 2 ) gen- 
erated using ideal avalanche photon detectors and the SSCS 
of amplitude a (solid curve) and the optimal fidelity with the 
corresponding regular SCS (dashed curve). The squeezing 
parameter of the initial squeezed state is r = —0.7. The pa- 
rameter of real transformation is T = 0.99. (b) The squeezing 
parameter of the target SSCS for which the fidelity is opti- 
mized. 



FIG. 6: Optimal fidelity F opt between the PSAS (a' a) gen- 
erated using ideal avalanche photon detectors and the SSCS 
of amplitude a (solid curve) and the optimal fidelity with 
the corresponding regular SCS (dashed curve). The squeez- 
ing parameter of the initial squeezed state is r — —0.7. (b) 
The squeezing parameter of the target SSCS for which the 
fidelity is optimized. The parameters of real transformations 
are T = 0.99, G = 1.01. 



B. Effects of imperfect detectors 

We have so far considered ideal avalanche photodetec- 
tors with unit quantum efficiency. However, in real ex- 
periments, detection efficiency is always limited. A real- 
istic detector with efficiency r\ can be modelled by a beam 
splitter with transmissivity r\ and vacuum at the idle port 
inserted in front of an ideal detector. The Wigner func- 
tion of the output state can then be represented as 



Wo Ut (xi,pi) 



1 



-^n [/(s: ~ , - 2/(sy (3 8 , 

-27(S;j + 47(E; v )] 

where the covariance matrice ££ ot incorporates the effect 
of the imperfect detection as 

s; ot = SE tot S + (l-»/)diag(0,0,i,i,i,i), 

E = diag(l, l.V^v^V^V 5 ?) ( 39 ) 
and the matrices E^, can be obtained from (l39l) and (l27l) . 



Finally, the fidelities can be arrived at using Eqs. (|36|) and 
([37]) . Figure [7] presents the optimal fidelity of the TPSS 
generated using avalanche photodetectors with rj = 0.6 
and a beam splitter of T = 0.99. Remarkably, under 
these realistic assumptions, the fidelity F w 0.95 can still 
be obtained. The fidelity for PSAS under the same con- 
siderations is plotted in Fig. [H] where the optimal fidelity 
drops down to F « 0.89. This confirms that two-photon 
subtraction is a better scheme to generate SSCSs of high 
fidelity, even though the amplitudes are smaller in this 



IV. SUCCESSIVE APPLICATIONS OF 
PHOTON SUBTRACTION 

We have shown that the SSCS is extremely well ap- 
proximated by the TPSS. We now show that this can 
be generalized to arbitrary TV-photon subtraction with 
N > 3. Namely, TV-photon-subtracted squeezed states 
(NPSSs) are good approximations of the SSCSs, which 
may be compared with the proposal by Fiurasek et al. to 
generate an arbitrary state by photon subtractions and 
displacements of a squeezed state [39| . We suppose that 
a beam splitter of transmmittivity T and an ideal pho- 
todetector is used to subtract N photons from a squeezed 
state. The Wigner function of the NPSS is then obtained 
as 38] 



B 2 N 

W N {(3) =M N exp ( - \p 2 r - ) J2 



^ o k\[(N-ky] 



where 



M 



N 



N/2 

E 



N-2k 



(/V-2fc)!(fc!) 2 ' 



fc=0 

x = (i-n)/(i+n), 

K = Tltanhrl 2 



(40) 



(41) 

(42) 
(43) 



and if ra [a;] is the Hermite polynomial. When N is odd 
(even), the NPSS should be compared with the odd 
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FIG. 7: (a) Optimal fidelity F op t between the TPSS (a 2 ) 
generated using realistic avalanche photon detectors and the 
SSCS of amplitude a (solid curve) and the optimal fidelity 
with the corresponding regular SCS (dashed curve). The 
squeezing parameter of the initial squeezed state is r = —0.7. 
The parameters of realistic transformations are T = 0.99 and 
r\ = 0.6. (b) The squeezing parameter of the target SSCS for 
which the fidelity is optimized. 



FIG. 8: (a) Optimal fidelity F opt between the PSAS (a'a) 
generated using realistic avalanche photon detectors and the 
SSCS of amplitude a (solid curve) and the optimal fidelity 
with the corresponding regular SCS (dashed curve). The 
squeezing parameter of the initial squeezed state is r — —0.7. 
(b) The squeezing parameter of the target SSCS for which the 
fidelity is optimized. The parameters of realistic transforma- 
tions are T = 0.99, G = 1.01 and r) = 0.6. 



(even) SCS. The fidelity between the NPSS and the SSCS 
can be obtained using Eq. (|35[) . It is nontrivial to obtain 
analytical expressions of the fidelity for an arbitrary N. 
As presented in Fig. [9l we numerically assess the fidelity 
F and plot it from N = 3 to N = 8 against the amplitude 
a of the SSCS. We suppose K is 0.6, i.e. r w 0.7 (6.1dB) 
when T = 0.99 and r = 0.75 (6.5dB) when T = 0.95. The 
degrees of squeezing, 6.1dB to 6.5dB, are experimentally 
achievable. Interestingly, the squeezing parameter of the 
SSCS is r' — 0.44 regardless of N. Our numerical cal- 
culation confirms this generalization up to W = 15. For 
example, the fidelity is as high as F > 0.999 for a = 2.52 
(a — 3) when N = 10 (N = 15) with the same squeezing 
parameter r' = 0.44 of the SSCS. 



V. REMARKS 

We have pointed out the twofold directions of the de- 
velopment for the generation of SCSs. The generation 
of SCSs with high fidelities (F > 0.99) and those with 
large amplitudes (a > 2) aims at practical applications 
of quantum information processing and macroscopic tests 
of quantum theory. The SSCSs, which were generated in 
recent experiments [30], may be a good alternative to 
the SCSs for the aforementioned purposes while their fi- 
delity in those experiments is yet to be improved. In this 
paper, we have studied how photon addition and sub- 
traction can be used to generate the SSCSs with high 
fidelities and large amplitudes. We have found that the 
single photon subtraction and subsequent addition with 
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FIG. 9: The optimal fidelity between the NPSSs with TZ(= 
Ttanhr) = 0.6 and the SSCSs of r' = 0.44 from N = 3 (left) 
to iV = 8 (right). The x axis represents the amplitude a of 
the target SSCS. The optimal fidelity is F > 0.999 regardless 
of the number of N. 



a squeezed vacuum state can cause the production of the 
approximate SSCS with F i=s 0.956 and a m 2 when the 
squeezing degree is about 6.1dB which is achievable using 
current technology. Furthermore, we show that N pho- 
ton subtraction may be used to generate the SSCSs with 
extremely high fidelities as F > 0.999. The amplitude of 
the SSCS increases as the number of the subtracted pho- 
tons gets larger. For example, a = 1.26 is obtained for 
two-photon subtraction while a = 1.68 for four-photon 
subtraction when the fidelity is F > 0.999 for both cases. 
We have assessed some experimental imperfections in 
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implementing photon addition and subtraction such as 
inefficiency of photodetectors and nonunit transmmittiv- 
ity of beam splitters. It has been shown that the fidelity 
of F w 0.95 can still be obtained with detection efficiency 
r\ = 0.6 and beam-splitter trasmmittivity T = 0.99 for 
two-photon subtraction to generate SSCSs. We believe 
that our work will immediately motivate experimental ef- 
forts to generate high-fidelity SCSs and large-amplitude 
ones that are useful for various quantum information ap- 
plications and fundamental studies. 
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